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Abstract 

We consider director configurations of cholesteric liquid crystal (CLC) cells with two 
plane confining substrates. Exact solutions of the Euler-Lagrange equations for out-of-plane 
orientations of the easy axes that correspond to inhomogeneous conical structures of CLC 
director are derived. We study dependence of the CLC twist wavenumber on the free twisting 
number assuming that anchoring energies at the substrates are either equal or different. In 
both cases this dependence is found to be generally discontinuous with hysteresis loops and 
bistability effects involved. For CLC cells with identical substrates and planar anchoring con- 
ditions the jump-like behaviour only disappears in the weak anchoring limit. Contrastingly, 
when the anchoring strengths are different, there is the finite value of anchoring below which 
the dependence becomes continuous. Another effect is the appearance of the gap between 
the adjacent twist wavenumber intervals representing locally stable director configurations. 
We calculate the critical value of anchoring asymmetry and present the results of numerical 
calculations. 
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1 Introduction 

In equilibrium cholesteric phase molecules of a liquid crystal (LC) align on average along a local 
unit director n(r) that rotates in a helical fashion about a uniform twist axis 0. This tendency 
of cholesteric liquid crystals (CLC) to form the helical twisting pattern is caused by the presence 
of anisotropic molecules with no mirror plane — the so-called chiral molecules (see for a recent 
review) . 
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The phenomenology of CLCs can be explained in terms of the Frank free energy 



Fb=^j d^x{Ki{V ■ n)2 + K2{n ■ V x n)^ 

+ [fi X (V X n)]2 + 2/i fi ■ V X n} , (1) 

where Ki, K2 and are the splay, twist and bend Frank elastic constants. As an immediate con- 
sequence of the broken mirror symmetry, the expression for the bulk free energy (|I]) contains the 
chiral term proportional to the chiral strength parameter h. This parameter gives the equilibrium 
value of CLC twist wavenumber, go = h/K2, which corresponds to the pitch Pq = 27r/go- In what 
follows the parameter qq will be referred to as the free twist wavenumber or the free twisting num- 
ber. So, if the twist axis coincides with the 2;-direction, the director field n = (cos qoz, sin qqz, 0) 
defines the equilibrium configuration. Since n and — n are equivalent in liquid crystals, periodicity 
of the spiral is given by the half-pitch, Po/2. 

Typically, the pitch Pq can vary from hundreds of nanometers to many microns or more, de- 
pending on the system. The macroscopic chiral parameter h (and thus the pitch) is determined 
by microscopic intermolecular torques |^ and depends on the molecular chirality of CLC con- 
sistuent mesogenes. In experiments it can be influenced either by variations in thermodynamic 
parameters such as temperature or by introducing impurities. 

So far we have seen that the twisting number of the ideal helical director configuration in 
unbounded CLCs equals the free twisting number go- This is no more the case in the presence 
of boundary surfaces or external fields that generally break the translational symmetry along the 
twisting axis. The case of planar CLC cell bounded by two parallel substrates exemplifies the 
simplest confining geometry and is of our primary concern in this paper. 

Director configurations in the CLC cell are strongly affected by the anchoring conditions at the 
substrates, so that, in general, the helical form of the director field will be distorted. Nevertheless, 
when the anchoring conditions are planar and out-of-plane deviations of the director can be 
neglected, the configurations still have the form of the helical structure which twist wavenumber, 
g, differs from go. Dependence of the twist wavenumber g on the free twisting number go has 
been previously studied in Refs. [|, |^, H and was found to be discontinuous. In Refs. [^, [l^ 



this jump-like behaviour was shown to manifest itself in abrupt changes of selective transmission 
spectra with temperature. Different mechanisms behind temperature variations of the pitch in 



CLC cells have been discussed in recent papers [0, Q 



Theoretical results of Refs. p, ^, ||, |T0|, |TT], |T2[ are related to the CLC cells placed between 



two identical substrates and the corresponding anchoring strengths were taken to be equal. In 
this paper we consider CLC cells with different anchoring energies at the confining surfaces. 
In particular, we find that sufficiently large asymmetry in anchoring strengths will render the 
dependence of the twist wavenumber on go continuous. 
The paper is organized as follows. 

Director configurations within CLC cells are analyzed in Sec. 0. Using the one-constant approx- 
imation we derive exact solutions of the Euler-Lagrange equations. These solutions are expressed 
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in terms of the Jacobian elhptic functions and correspond to inhomogeneous conical orietational 
structures. The case of the planar anchoring conditions is considered in Sec. ^ where we present 
both analytical and numerical results about the dependence of the twist wavenumber on the free 
twisting number. Concluding remarks are given in Sec. ^. 



2 Orientational structures in CLC cell 

We consider director configurations in a CLC cell of thickness 21 sandwiched between two parallel 
plane substrates located at 2; = — / and z = I, so that the 2;-axis is normal to the confining surfaces. 
We shall also assume that the symmetry with respect to in-plane translations is unbroken and the 
director field describing orientational structures does not depend on the coordinates x and y. Its 
parameterization in terms of the angles and 6 is taken in the form: 

h{z) = {cos(j){z) cos 9 (z), sin (f){z) cos 9 (z), sin 9 (z)) . (2) 

Substituting the parameterization into Eq. (|1|) and using the one-constant approximation, 
Ki = K, we have the following expression for the bulk free energy per unit area: 



K 
2 



j [{9'f + {(!)' cos^ 9- q^f + {(t)' sin9 cos 9f]dz, (3) 



where the prime stands for derivative with respect to z. 

In order to specify the anchoring conditions at the boundary surfaces 2; = ±/ in terms of 
the anchoring strengths W± and the vectors of easy orientation e± we shall write the surface 
contribution to the free energy as an anchoring energy taken in the form of Rapini-Papoular 
potential: 



l-(n(/)-e+)^ +^ l-(n(-/)-e_)^ . (4) 



The phenomenological anchoring parameters W± characterize the strength of interaction between 
CLC molecules and the surfaces, while the preferred orientation of the molecules at the substrates 
is defined by the easy axes e±. 



2.1 Solutions of the Euler-Lagrange equations 

The free energy functional (|^) only depends on the derivative 0', so that the corresponding Euler- 
Lagrange equation can be integrated to yield the following equation: 

0' = go(l + ccos-'0), (5) 

where c is the integration constant. Variation of the functional with respect to 9 will give the 
second Euler-Lagrange equation that can be simplified by using Eq. (|^). The result is as follows 

9" = -qlsin9cos9{c^cos'H-l). (6) 
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Figure 1: Real- valued solutions of Eq. (^) exist only if the integration constants c and d define a 
point in the c-d plane that is above the curves: ci = + 1 at |c| > 1 (solid line) and d = 2\c\ aX 
|c| < 1 (dashed line). 



The first integral of Eq. (P) is not difficult to find and is given by 

{e'f + g2(c2 cos-2 e + cos^ e) = qld , (7) 

where d is the integration constant. From Eqs. (|^) and (P) we have the bulk free energy in the 
form ^ 

2F^,/K = ql j {d-2cos^e + l)dz. (8) 

We can now proceed with solving Eq. (|^) that can be conveniently rewritten in the form: 

[(sine)f = -go'(cos2 9 - c_)(cos2 9-0+), (9) 

where 2c± = ± \/ — Ac^. This equation will have solutions only if the values of c_ are non- 
negative and do not exceed the unity, < c_ < 1. This immediately places the restriction on the 
parameter d: d > max(2|c|, + 1). Fig. |l] shows this region in the c-d parameter plane. 
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Figure 2: Dependencies of /3 (= 2ql — Acp) on /5o (= 2go^ — A0) at = = W for various 
values of the dimensionless anchoring energy parameter w (= Wl/ K): (a) strong anchoring hmit, 
w oo] {h) w = tt; (c) w = n/2; (d) w = 0.25. 



In the upper part of the region, d > + 1, we have c+ > 1 and solutions of Eq. can be 
expressed in terms of the Jacobian elliptic functions [|1^ as follows 

sin 6* = a cn(6go2; + ^o|"^ "'^) ; (10) 

where = 1 — c_, 6^ = c+ — c_, m = and 6q is the integration constant. Below the parabola 
d = + 1, where 2\c\ < d < + 1 and c+ < 1, these solutions are given by 

sin 6* = a dn(ago-z + 6*01^) . (11) 



In addition, there are homogeneous solutions that are independent of z: sin 6* = — c_ and 
sin^ = ±Y^1 — c+ (at c+ < 1). These solutions correspond to a conical helix in which the director 
makes an oblique angle with respect to the helix axis. There are no spatially varying solutions at 



5 




Figure 3: The curves representing the plot of the function 7+(/3) and 7-(/3) are shown as thick 
sohd and dashed hues, respectively. The points located at the intersection of the curves and the 
horizontal straight line 7 = /^o give the roots of Eq. (^). The value of /3o is (10 + l/2)7r (thin solid 
line) and is (10 + l/2)7r±w (thin dotted lines). Two cases are illustrated: (a) w+ = w_ = w = 7r/2; 
(b) w. = 10.0, w+ = 2.19, w = tt/2 (see Eq. (g^)). 



the boundary of the region depicted in Fig. |1]. In this case the pretilt angle does not depend on z 
and we have 

. n Jo, 6 = c2 + l and |c| > 1, 

smfc' = < , (12) 

[±VT^, 6 = 2|c| and |c| < 1. ^ ^ 

Given the pretilt angle as a function of z, integrating of Eq. (^ will give the azimuth angle 0: 

<P{z) = qoZ + cJ^ ^3-^ + <^o, (13) 

where 0o is the twisting angle in the middle of CLC cell. 

In order to compute the integration constants c, d, 6q and 0o "we need to know the boundary 
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Figure 4: The free energy /+ (solid line) and /_ (dashed line) of stable configurations with the 
half-turn numbers between 10 and 14 as a function of f3o computed from Eq. (0) by using Eq. (|23|) 
for various values of the anchoring energy parameters (a) w = tt/2; (b) w = 7C. 



conditions. The standard variational procedure will yield the following result: 

Kqoc = ±W±cos9{n- e±){h^ ■ e±)\,=±i, (14) 
KO'U^^i = ±W±ih ■ e±)(n, ■ e±)|,=±, , (15) 

where n,^ = (— sin 0, cos0, 0) and = (— sin 9 cos0, — sin 6* sin 0, cos 6'). 



3 Pitch wavenumber versus free twisting number 

In this section we consider the simplest case of an ordinary spiral director configuration in which 
the pretilt angle equals zero, sin 6' = and d = + 1. This case occurs when the anchoring 
conditions are planar and both of the easy axes are parallel to the confining surfaces. We direct 
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Figure 5: Dependencies of (3 on /3o calculated at W- = 10.0 for various values of the parameters a 
(= 2w-W+/{w- — and w (see Eq. (a) w = vr (a = 11.8 and w+ = 3.705); (b) w = 11/2 

{a = 5.6 and w+ = 2.19); (c) a = 1.0 («;+ =Wc = 0.476); (d) a = 0.5 {w+ = 0.25). 



the X-axis along the vector e_, while the easy axis at the upper substrate is rotated through the 
angle A0: e+ = (cos A(/),sin 0). 
From Eq. (0) we have 

^ = qz + (f)Q, (16) 

where q = qo{l + c) is the twisting wavenumber. The boundary conditions (|T^) can be conveniently 
rewritten in the form: 

p-po = -w± sin(/3 ± a) , w± = W±l/K , (17) 
p = 2ql- A0, a = 200 - A0, po = 2qol - A0 . (18) 

Eqs. ([T7| ) define stationary points of the free energy written as a function of the parameters a and 
4/F(a, /3) /K ={13- - w+ cos(/? + a) - w_ cos(/3 - a) + (w_ + w+). (19) 
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Figure 6: The free energy (solid line) and /_ (dashed line) of the configurations with the 
half-turn numbers between 10 and 14 as a function of (3 computed from Eq. (p^). Thin lines 
represent the energy of unstable configurations. Two cases are shown: (a) w_ = Wj^ = w = 7r/2; 
(b) w_ = 10.0, w+ = 2.19 (d = 5.6 and w = tt/2). 



There are two additional conditions for a stationary point to be a local minimum of the free energy 
surface (|l^). These conditions ensure local stability of the corresponding director configurations 
and can be written in the following form: 

A = w+ cos(/? + a) +w_ cos(/? - a) > , (20) 
H = A + 2w+w_ cos(/3 + a) cos{(3 - a) > 0. (21) 

From Eq. (0) we can now relate the angles a and (3 through the following equation 

a = arctan[etan/3] + irk , e = — , (22) 

where the number k is integer, /c G Z, that is the number of half-turns of the spiral. Substituting 
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the relation ( p2D into Eq. ( p!?]) will give a pair of transcendetial equations for (3: 

Po = 7±(/3) = P ± 1(7+ sin (/5 + arctan[etan/3]) . (23) 

The right hand side of Eq. (|23|) is 7+(/5) (7_(/?)) when the number k is even (odd). Solving Eqs. (|23|) 
will provide, in general, various values of /?. The values that meet the stability conditions (PI1| ) 
and (|2T|) correspond to metastable configurations of CLC. The equilibrium director structure is 
determined by the solutions of Eq. ( PBj ) that give the least value of the free energy. From Eqs. (|T^) 
and (^3]) the free energy can be expressed in terms of the angle f3 and is given by 

4:IF±{(3)/K = f±{(3) = wlsin'^ =f {w+costlj+ + w-cosifj-) + (w_ + w+), (24) 
V'i = /3 ± arctan[etan/5] . (25) 

Similarly, we can use Eq. (|22|) to rewrite the stability conditions as follows 

A± = ±(u'+ cos 4'+ + W- cos > , (26) 
H± = A± + 2w^w^ cosip+ cosip- > 0. (27) 
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It is now our task to study the dependence of the pitch wavenumber q on the free twist 
wavenumber go? which is proportional the chiral parameter, go = h/K, at different anchoring 
conditions. Equivalently, we concentrate on the dependence of (3 on (3q (see Eqs. (|T^) and 
This dependence can be thought of as a sort of dispersion relation. 



3.1 Strong anchoring limit 

First we consider the case of the strong anchoring limit, W± = W oo^ that can be readily 
treated without recourse to numerical analysis. In this case the boundary conditions require the 
director at the substrates to be parallel to the corresponding easy axes, n(±0 || e±. This imposes 
the restriction on the values of g, so that g takes the values from a discrete set. This set represents 
locally stable director configurations labelled by the half-turn number m as follows 

2ql = A0 + 7rm, m G Z . (28) 
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Substituting the values of (3 from the relation (pSf ) into the first term on the right hand side of 
Eq. ( plf ) will define the equilibrium value of m as the integer that minimizes the distance between 
Tim and 2go/ — A0. The resulting step-like dependence is depicted in Fig. 



3.2 Equal anchoring energies 

When the anchoring strengths at both substrates are equal, = = W, the right hand side 
of Eq. (^) is P ± w sin P {w = Wl/K) and = V'- = /5- So, as is illustrated in Fig. 0(a), we 
need to find the intersection points of the horizontal line 7 = /3o and the curves 7 = 7±(/3) in the 
/? — 7 plane. The stability conditions > and > {A_ > and H_ > 0) require the 
values of f3 to be ranged between (m — l/2)7r and (m + l/2)7r, where m is even (odd) integer equal 
to the half-turn number. The function 7+ (7^) monotonically increases on these intervals, so that 
the value of (3q runs from (m — 1/2)7? — w to {m + l/2)7r + w on the interval with the half-turn 
number m: [{m — l/2)7r, (m + l/2)7r]. For each number of m we have the monotonically increasing 
branch of the (3 vs (3o dependence. The branches with m ranged from 10 to 14 for different values 
of the dimensionless anchoring energy parameter w are depicted in Fig. ^ It is illustrated that 
the /?o-dependence of (3 will be discontinuous provided the anchoring energy is not equal to zero. 
Fig. ^ shows that the jumps tend to disappear in the limit of weak anchoring, where the anchoring 
energy approaches zero, w —>■ 0. 

Similar to the case of strong anchoring, at (3o = vr/2 + vrm we have two different roots of 
Eq. (p3D with k = m and k = m + 1 that are equally distant from Pq and are of equal energy. 



According to Refs. |I0|, |Tl], [12| it can be assumed that the jumps actually occur at the end 
points Pq = tt/2 + nk ±w, where the configuration with m = + 1/2 =]= 1/2 becomes marginally 
stable {A± = H± = 0) and loses its stability. In other words, the system needs to penetrate the 
barrier separating the states with different half-turn numbers. As it is seen from Fig. 0, in this 
case the upward and backward transitions with Am = ±1: k k + 1 and k + 1 —>■ k occur at 
different values of jSo: jS^'^ = 7t/2 -\- Tvk + w and l3^^ = n /2 + nk ~ w, correspondingly. It means 
that there are hysteresis loops in the response of CLC to the change in the free twisting number. 

When the anchoring energy is small and w < 7r/2, there are only two configurations at the 
critical point Po = the marginally stable configuration with m = k and the equilibrium one 
with m = k + 1. In this case Eq. (|2^) has at most two roots and the jumps will occur as transitions 
between the states which half-turn numbers differ by the unity, |Am| = 1. 

Fig. ^ illustrates that for w = 7r/2 we have two marginally stable structures of equal energy 
with m = k and m = k + 2. Being metastable ai n/2 < w < tt this newly formed structure 
and the configuration with m = k + 1 will be degenerate in energy at w = vr. So, as is shown in 
Figs. 03 and ^d, both transitions k k + 1 and k ^ k + 2 are equiprobable, so that we have 
the bistability effect at the critical point under 10 = 11. For it < w < 2tt the configuration with 

(k) 

m = k + 2 gives the equilibrium director structure at Po = P\ . 

It is not difficult to see that further increase of the anchoring energy involves the configuration 
with m = k ±1 into the transition when the parameter w is passing through the value (/ — 3/2)7r. 
The jump to the equilibrium state will require the half-turn number changed from k to k + Am, 
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where | Am| = / at {I — 1)7t < w < Itc . There are a number of transient metastable configurations 
involved in transitions with I Am I < /. 



3.3 Different anchoring energies 

When the anchoring energies at the surfaces are different, W- ^ W+, sinip^ on the right hand 
side of Eq. ( p3D equals zero at P = tt/2 + nk and, as is seen from Fig. ^d, we have additional 



intersection points of the curves 7+(/3) and 7-(/3). The stationary points, where the derivative of 
7± with respect to f5 equals zero, represent the local maxima and minima of 7-1- and are located 
at = 7r/2 + vrfc ± A/3. The value of A/5 can be calculated by solving the following equation 

/.^ r . ^nx 1 + (e^ - 1) COS^ A/? , , 

+ e) sin(A/3 - e arctan[e cot A/3]) = - ^ _ , (29) 

where A(3 G [0,7r/2]. 

From the stability conditions > > 0) the values of jS for stable configurations 
fall between the stationary points (m — l/2)7r + A/5 and (m + l/2)7r — A/3, where the half-turn 
number m is even (odd) integer. The function 7+ (7.) monotonically increases and (3q varies from 
(m — 1/2)7!- — w to (m + 1/2)tt + w on these intervals. The dimensionless parameter w, as opposed 
to the case of equal anchoring energies, is now given by 



w 



w+cos(A/5 - earctan[ecot A/5]) - A/5. (30) 



Clearly, we can now follow the line of reasoning presented in the previous section to find out 
the results concerning hysteresis loops and bistability effects that are similar to the case of equal 
anchoring strengths. There are, however, two important differences related to Eqs. (^9]) and (^Ol). 

If A/5 7^ 0, the intervals of /3 representing stable director configurations are separated by the 
gap of the length 2A/5. The presence of this gap is illustrated in Figs. ^ and |]b. Fig. shows 
the gap between stable branches of the dependence of the free energy on /5. The values of (3 that 
are within the gap represent unstable configurations ("forbidden" states of CLC). The graph of 
the A/5 vs w+ dependence is presented in Fig. |^ and indicates that the gap disappears in the limit 
of equal energies, w+ = w_. In addition, referring to Fig. |^, there is a critical value of w+ below 
which A/3 equals zero. 

It can be shown that Eq. (^) has the only root, /5 = /5o at /5o = 7r/2 + 7rA;, under the anchoring 
energies meet the following condition: 

a^^^^<l. (31) 

In this case the gap disappears and the dependence of /5 on Pq becomes continuous in the manner 
indicated in Figs. ^ and ^i. Given the value of ty_ the relation (^) yields the threshold value 
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for the anchoring strength at the upper substrate. The critical anchoring energy parameter Wc is 
given by 

Wc = — — — . (32) 

As illustrated in Fig. the dependence of w on w+ is approximately linear for large values of 
W- and w goes to zero at the critical point = Wc- 



4 Conclusion 

In this paper we have studied how the pitch wavenumber of CLC cell depends on the free twisting 
number for the cases in which the anchoring strengths at the substrates are either equal or different. 
It is found that in both cases this dependence is generally discontinuous and is characterized by 
the presence of hysteresis and bistability. But the difference in the anchoring energies introduces 
the following two effects: 

(a) the jump-like behaviour of the twist wavenumber is suppressed under the anchoring strength 
at one of the substrates is below its critical value; 

(b) the twist wavenumber intervals of locally stable configurations with adjacent numbers of the 
helix half-turns are separated by the gap where the structures are unstable. 



The part of our analysis, presented in Sec. 3.2, relies on the assumption that the transition 



between configurations with different half-turn numbers occurs when the initial structure loses 
its stability, so that its pitch wavenumber no more corresponds to a local minimum of the free 
energy surface. The result is that the stronger the anchoring, the larger the change of the half- 
turn number (and of the twist wavenumber) needed to reach the equilibrium state. So, whichever 
mechanism of relaxation is assumed, the metastable states certainly play an important part in the 
problem. Dynamics of the transitions is well beyond the scope of this paper and it still remains 
a challenge to develop a tractable theory that accounts for director fluctuations, hydrodynamic 
modes and defect formation. 

Analytical results of Sec. |] can be used to study the effect of pretilt angles at the confining 
surfaces. This problem requires a more complicated analysis that will be published elsewhere. 
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